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ABSTRACT 

Although the driving mechanism a cting in (3 Cephei pulsators is well known 
(e.g. iDziembowski fc Pamvatnvkhl fl993) . problems concerning identification of am- 
plitude limitation mechanism and non-uniform filling of the theoretical instability 
strip, remain to be solved. In the present analysis, these problems are addressed by 
non- linear modelling of radial pulsations of these stars. In this approach radial modes 
are treated as representative for all acoustic oscillations. 

Sev eral models o f differ ent masses and metallicities were converged to limit cycles 
through IStellingwer f (1974) relaxation technique. Resulting peak-to-peak amplitudes 
are of order of AV = 0.3 mag. Such amplitudes are significantly larger than those 
observed in (3 Cephei pulsators. Assuming that all acoustic modes are similar, we show 
that collective saturation of the driving mechanism by several acoustic modes can easily 
lower predicted saturation amplitudes to the observed level. Our calculations predict 
significant decrease of saturation amplitudes as we go to high mass/high luminosity 
models. However, this effect is not strong enough to explain scarcity of high mass 
(3 Cephei variables. A possible weakness of the collective saturation scenario is that 
estimated line-broadening, resulting from excitation of many high-/ modes, might be 
higher than observed in some of the (3 Cephei stars. We argue that this difficulty can 
be overcome by allowing g-modes to participate in the saturation process. 

We also discuss robust double-mode behaviour, encountered in our radiative mod- 
els. On a single evolutionary track we identify two double-mode domains with two 
different mechanisms resposible for double-mode behaviour. The non-resonant double- 
mode domain separates first overtone and fundamental mode pulsation domains. The 
resonant double-mode domain appears in the middle of the first overtone pulsation 
domain. Its origin can be traced to the 2uj\ = 0^0 + ^2 parametric resonance, which 
destabilizes the first overtone limit cycle. 

Key words: hydrodynamics - stars: early type - stars: oscillations - stars: variables: 
other. 



1 INTRODUCTION 

13 Cephei stars are main sequence early B-type pulsators. 
Their radial and non-radial p-modes and low-order g- 
modes are driven through the c lassical /t-mechanism act- 
ing in the metal opacity bump llMoskalik fc Dziembowski 



1992 ; IDziembowski fc Pamvatnvkh 19931 : Gautschv fc Saic 
19931 ). Most of /3 Cephei stars are multiperiodic with typ- 
ical periods ranging from 3 to 7 hours, and peak-to-peak 
l/-band amplitudes usually below 60mmag. For observa- 
tio nal review of astrophysical pro perti es of j3 Cephei stars 



Sterken fc Jerzvkiewiczl (|l993l ) and IStankov fc Handled 

The latter paper contains the catalogue of 93 con- 
firmed galactic f3 Cephei pulsators. This list should be sup- 



plemented by 1 9 new objects found in ASAS-3 survey data 
jPigulskill2005l : iHandlerl feoOS). Most of these newly discov- 
ered objects have high amplitudes in comparison with pre- 
viously known variables (cf. Table [T]). This is a selection 
effect due to automatic algorithm used in the ASAS sur- 
vey to extract variable objects from several million stars be- 
ing monitored. Only stars showing an excessive dispersion 
are s ubject do detailed periodogram analysis (|Poimanskil 
2002). Thus, low amplitude pulsators are not present in the 
ASAS catalogues. The ongoing periodogram analysis of all 
the ASAS-3 stars of proper spectral types at least triple s 
the number of detected (5 Cephei pulsators l|Pigulskill2007l ). 
Several f3 Cephei variables were also found in the Magellanic 
Cloud s ( Pieulski fc Kolaczkow ski 2002; Kol aczkowski et al.l 
l2004bl . [200rf )~ 
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Although the driving mechanism acting in f3 Cephei 
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stars is well known, several problems remain to be solved. 
These concern observational boundaries of f3 Cephei in- 
stability strip and question of amplitude limitation mech- 
anism. The theoretically predic ted instability strip is not 
completely filled with stars (e.g. IPamvathvkhl 1 1993 . Fig. Q] 
this paper). Particularly, linear analysis predicts unstable 
p-modes in slightly ev olved massive O-type stars. De spite 
of observational effort l|Pigulski fc KoIaczkowskilll998l ). un- 
til very recently we have not known any firm pulsat or of this 
kind. Three candidates were found only recently l|Pigulskil 
2007). Also, lack of stars at low-mass end of the theoret- 
ical instability strip remains to be explained. Considering 
detected frequencies, about 38 per cent of galactic (3 Cephei 
pulsators are monoperiodic. For the purpose of further com- 
parison, in Table 1 we have listed the monoperiodic vari- 
ables for wich V-band amplitudes are available. Selectio n 



of mo noperiodic variables was based on lStankov fc Handler 
(2005 ) catalog and p ublications on ASAS variables IjPigulski 
120051 ; iHand lcr 2005). Peak-to-peak amplitudes of these vari- 
ables are taken from various publications (references are 
given below the Table). The Table contains both field and 
cluster variables and is based on analysis of observational 
data of differing length and quality. Obviously, monoperiod- 
icity of variables listed in Table 1 is a statement reflecting 
our current knowledge. 

Although significant fraction of /3 Cephei variables ap- 
pears to be monoperiodic, it is very likely, that almost 
all of the P Cephei stars are in fact multiperiodic, with 
amplit udes of additional mode s below current detection 
limit (|Stankov fc Handled 120051 ). This hypothesis is sup- 
ported by intensive observation campaigns of individual 
variables, leading to detection of mor e pulsation modes 
with decreasing detection threshold (e.g. Ijerzvkiewicz et aU 
120051 ). We also note, that t wo st ars claimed as monoperi- 
odic in IStankov fc Handlei] (120051 ) catal ogue, are now con- 
firmed multi mode pulsator s (8 Cet - lAerts et al.l 120061 ; 



HD 164340 - Pigulski 2005) . From observ a tional point of 



view Ibaszvhs ka-Daszkiewicz fc Niemczural l|2005h showed 
that multiperiodic variables have a larger metal abundance 
in the photosphere than monoperiodic ones. Linear calcu- 
lations also show that for higher metal abundance we g et 
more unstable modes (|Dziembowski fc Pamvatnvkhll 19931) . 

Non-linear effects are responsible for the amplitude lim- 
itation and thus, non-linear calculations may shed light on 
problems concerning observed instability strip and insta- 
bility saturation. So far full non-linear calculations can be 
done for radial modes only. R adial modes are not r are in 
the observed f3 Cephei stars (|Aerts fc De Catl 120031 ). Ex- 
cept of BW Vulpeculae (which is an extreme case), radial 
modes have similar amplitudes to non-radial modes. Thus, 
we assume that a radial mode is a typical representative of 
all acoustic modes. This approximation is applicable to p- 
modes only and not to g-modes as their physicall nature and 
thus saturation properties are different. 

First non-linear models of (3 Cephei va r iables 
were calculated for B W Vul JPesnell fc Cox! Il980l ; 
iMoskalik fc Buchledll99l ; iFokin et al.ll2004l ). This monope- 
riodic pulsator exhibits most extreme light variability 
among known (5 Cephei stars, with peak-to-peak a mpli- 
tude approximately 0.2 mag in V. Moskali k~fc~Buchlen con- 
cluded that BW Vul is radial, fund amental mo d e pu l- 
sator, which was further confirmed by lAerts et alj l|l995h . 



IMoskalik fc Buchleij were able to reproduce the observed 
amplitude, which leads to conclusion that in this case insta- 
bility may be saturated by a single mode. However, BW Vul 
is an extreme case, since most /3 Cephei stars have much 
smaller amplitudes. F or example 7 Peg is a monoperiodic 
radial mode pulsator l|Aerts fc De Catll2003f ) with F-band 
amplitude of 17mmag (Table [TJ. In the present analysis, 
conclusions about instability saturation are drawn through 
comparison of calculated non-linear limit cycle amplitudes of 
P Cephei models with observed amplitudes of /3 Cephei stars. 
Our preliminary results were presented at the Vienna Work- 
shop on the Future of A stero seismology l|Smolec fc Mo skalik 
l2007f) . 

This paper is organised as follows. In Section 2 we 
present our hydrocodes and model construction and anal- 
ysis. In Section 3 we discuss sensitivity of calculated ampli- 
tudes to several numerical parameters (artificial viscosity, 
numerical mesh) . Section 4 contains main results of this pa- 
per concerning mode saturation. In Section 5 interesting case 
of double-mode pulsation in non-convective stellar models is 
discussed. Conclusions are summarized in Section 6. 



2 CONSTRUCTION OF MODEL SEQUENCES 

Evolutionary main-sequence models were computed using 
Warsaw-New Jerse y evolutionary code (for details see e.g. 
IPamvathvkhl [l999l . and references therein). Rotation and 
overshooting from the convective core were disregarded. As 
will be justified in Section 4.4, convection in model en- 
velopes may be neglected since /3 Cephei stars are too hot 
to have an effective energy transfer by convection in stel- 
lar envelope. Several tracks with metallicities Z = 0.015 
and Z — 0.02, and masses 7-14 Mq were calculated. Ad- 
ditional high-mass track (20 Mq, Z — 0.02) was also com- 
puted. The OPAL opacities were used in all evolutionary 
calculations. Computed evolutionary tracks are presented 
in Figure [T] For comparison we also plot galactic /3 Cephei 
variables (circles in Fig. [TJ. Their bolometric luminosities 
and effective temperatures were obtained using Stromgren 
photometry indices collected in IStankov fc Handlerl (|2005l ) 
catalogue. E ffective temperatures were obtained using cal- 
ibrations of iMoon fc Dworetskvl l|l985l ) (uvbybeta code 
written by Moon and modified by Napiwotzki). My and 
bolometric corrections wer e derived via calibrat ions of 
Balonafc Shobbrookl (|l984l l and from the work of iFloweil 
(1 19961 ). respectively. Monoperiodic variables are plotted as 
open circles. They do not form any distinct group. 

Bolometric luminosities and effective temperatures, re- 
sulting from evolutionary calculations, were used in con- 
struction of model envelopes, that were subject to further 
radial pulsation analysis. Static envelopes were calculated 
using 120 zone mesh, with zone mass increasing geometri- 
cally inward. The temperature of the 70th zone from the 
surface was fixed to 2.5 x 10 5 K. The zone mass ratio above 
this anchor was the same for all models, while below the an- 
chor it was iterated to match the desired bottom boundary 
temperature, Tin, equal 3 x 10 7 K. Described mesh assures 
reasonable resolution in a shock region developing in exterior 
zones of non-linear models. 

All models were subject to linear pulsation analysis 
and selected models were subject to non-linear calcula- 
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Table 1. Monoperiodic j3 Cephei variables, for which l/-light am- 
plitudes are available. Last seven stars are new variables detected 
in ASAS-3 catalogue. Amplitudes are peak-to-peak values. 



Star 


period [d] 




Ref. 


7 Peg 


ri 
u 


101O 


0.017 


1 


NGC 663 v4 





1940 


0.040 


2 


NGC 869 v692 





1717 


0.019 


3 


NGC 869 v992 





1326 


0.003 


3 


NOC 884 v2299 





3190 


0.016 


4 


4 CMa 





2096 


0.034 


I 


KK Vel 





2157 


0.050 


5 


NOC 4755 v301 





1305 


0.010 


5 


NHP 47^ v9D1 





1825 


0.010 


5 


nhp 47^^ v9i n 

I^VjVj 'ilUO V^J_U 





0933 


0.007 


5 


Lup 





1774 


U.Uz < 


i 


8 Lup 





1655 


0.0035 


7 


vl449 Aql 





1822 


0.074 


8 


NGC 6910 vl4 





1904 


0.017 


9 


NGC 6910 v24 





1430 


0.009 


9 


BW Vul 





2010 


0.2 


1 


HN Aqr 





1523 


0.032 


10 


CPD-62°2707 





1417 


0.110 


8 


HD 133823 





1760 


0.106 


8 


CPD-46°8213 





2242 


0.129 


8 


HD 328906 





1776 


0.084 


8 


HD 152477 





2650 


0.071 


8 


BD — 14° 5047 





2402 


0.088 


8 


HD 100495 





1685 


0.037 


11 



Periods are taken from [ S tankov fc Handler! ll2005h. Ampli- 



2001 


; (3) Krzcsinski et al. 


19991; (4) iKrzesihski & Pigulskil 19971; 


(5) 


Cousins! 11982 


; (6) 


Stankov et al.l |2002|; (7) Ishobbrookl 


1972 


; f8) Pieulsk 


2005; (9) Kolaczkowski et al.l l2004al; (10) 


Kilkenny & van Wyk 1990; (11) Handler 2005. 



Notes on individual stars: 8 Lup - two periods are claimed in the 
literature, howeve r photometric period is li kely an alias of the 
spectroscopic one dStankov fc Handler! I2OO5I) ; KK Vel - two fre- 
quencies / and 2/ are detected in this variable. However, it is not 
clear if 2/ is the harmonic frequenc y or an independen t, nonlin- 
early coupled mode as suggested by lAerts et al.l dl994T) . Even in 
the latter case, KK Vel can be treated as a monoperiodic variable, 
similarly to the bump Cephcids. 



tions. We used radiative Lagrangean hydrocodes, suitable 
for radi al pulsation anlysi s. These codes are essentially 
th ose of IStellingwerJ (Il975l) . with modifications described 
in iKovacs fc Buchlen |l98ct) . In the linear non-adiabatic 
(L NA) code we used canned eigenvalue solver as suggested 
bv lGlasner fc Buchlerl l|l995t ). We assumed ~ exp(crf) depen- 
dence of the perturbed quantities, where a is the complex 
eigenvalue. We define the linear growth rates of the modes 
as fractional growth of the kinetic energy per pulsation pe- 
riod: 7 = 4n?R.(a) / ui , where: ui = Q(<j), is mode frequency. 
Analytical equation of state (Stellingwcrf 1982), and most 
recent version of OPAL opacities l lglesias fc Rogers! Il996l ) 
were u sed. Opacity table s were generated for solar mix- 
ture of iGrevesse fc Noels! (|l993h . For some s equences OP 
opaci ties $Scaton 20051) with Seaton mixture l|Seaton et al.l 
1 19941 ) were used. Linear stability results are presented in 
Figure [T] Fundamental and first overtone instability strips 
(70 > 0, 71 > 0) are enclosed by dashed and dotted lines, 
respectively. First overtone is unstable for Z — 0.02 only. 



Higher-order overtones are always stable. Observed instabil- 
ity strip runs approximately parallel to the Zero- Age Main 
Sequence (ZAMS) and the Terminal-Age Main Sequence 
(TAMS) lines. Theoretical instability strip extends to post 
main sequence phase, but since evolution becomes very fast 
shortly after the TAMS, the TAMS is often considered as 
the effective red edge of the pulsation instability domain. 

For selected models, full amplitude non- linear limit cy- 
cles (periodic finite-amplitude pulsations) of fundamental 
and first overtone modes were calcula t ed using relaxatio n 
scheme dBaker fc von Sengbuschl Il969l ; IStellingwerJ I1974T I. 
The relaxation scheme provides also inform ation abou t 
limit cycle stability (Floquet analysis, see e.g. lined Il956l ). 
The same numerical mesh, difference scheme and opac- 
ity/equation of state subroutines were used both in LNA 
and in hydrocodes. Shocks were treate d with an artificial 
viscosity formula of IStellingwerJ |l975j), with artificial vis- 
cosity parameters Cq = 4, and cutoff parameter a = 0.06. 
The choice of viscosity parameters will be discussed in Sec- 
tion 3.1. In non-linear iterations we used 400 time steps per 
pulsation period and in all cases it was enough to obtain 
reliable amplitude. However, for evolved fundamental mode 
models with shocks, up to 2000 time steps per period were 
necessary to obtain reliable values of Floquet stability coeffi- 
cients. Since most of the observational data are in a V-band, 
theoretical lig ht curves were transformed to F-band through 
iKuruczl (|2006j ) static atmosphere models. This was done by 
applying a bolometric correction at each pulsation phase, 
however omitting the shock phase, where this quasi-static 
approximation is not applicable. 



3 SENSITIVITY OF CALCULATED 
AMPLITUDES 

Since conclusions of this paper rely on calculated limit cy- 
cle amplitudes, it is necessary to check how reliable they 
are. It is well known from classical pulsators modelling, that 
for a given model (M, T B s, L, X, Z), limit cycle satura- 
tion amplitudes may vary dependin g on artificial viscosity 
and numerical mesh (|Kovac s 1990). In this Section these 



effects are shortly analysed. All calculations are performed 
for bolometric light curves, since they are more smooth than 
V-band curves and void of transformation uncertainties (as 
we found, simple linear relation of type Av = a ■ Abo\ exist, 
thus results concerning amplitude uncertainty are the same 
for V-band curves). Tests were done for the fundamental 
mode pulsations of three models, denoted by A, B and C in 
Fig.[j] All models have M = 11 M Q , X = 0.7 and Z = 0.02. 
Models A and B are close to the boundary of the instability 
strip, while model C is an evolved model approaching the 
TAMS. In all calculations 400 time steps per pulsation pe- 
riod were used. Test calculations for other models and for 
the first overtone mode showed, that results presented below 
are representative. 



3.1 Sensitivity to artificial viscosity 

In case of radiative models of classical pulsators, artificial 
viscosity parameters determine the limit cycle amplitude, 
and must be chosen to satisfy the observational constraints. 
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Figure 1. Evolutionary tracks and LNA stability results for (3 Cephei models with Z = 0.02 (left panel) and Z = 0.015 (right panel). 
Dashed and dotted lines enclose the fundamental and first overtone instability strips, respectively. Circles represent observational data: 
open circles — monoperiodic variables, filled circles — multiperiodic variables. Diamonds refer to positions of test models, A, B and C, 
described in Section 3, while crosses refer to convective models described in Section 4.4. 



Fortunately, this is not the case for /3 Cephei stars. Calcula- 
tions performed for test model C show, that for sufficiently 
high cutoff parameter (a > 0.04) unphysical interior dissi- 
pation, which affects the amplitude, is eliminated. Artificial 
viscosity turns on only in the exterior shock region, with 
almost no effect on limiting amplitude. This is illustrated 
through non-linear work integrals (Fig. [2]) and bolometric 
light curves (Fig. [3} obtained with different cutoff parame- 
ters. With q = 0.02 (left panel of Fig.[5| one note significant 
negative contribution of artificial viscosity work in radiative 
damping zones (around zone 40), which leads to decrease of 
amplitude (Fig. [3} • This dissipation is unphysical since no 
shock is generated here. Thus, we need to increase the cutoff 
parameter in order to eliminate this unphysical dissipation. 
For q = 0.04 (right panel of Fig. [2Jl artificial viscosity con- 
tributes only in the exterior, shock zones. As it is well visible 
in Fig. [3] with increasing a, amplitudes converge to constant 
value. Indeed, light curves obtained for a £ (0.04, 0.1) over- 
lap and for clarity are not presented in Fig. [3] Also keep- 
ing a — 0.06 fixed and varying Cq parameter in a range, 
Cq £ {2, 3, 4, 5, 6}, resulting light curves are indistin- 
guishable. This justifies our choice: (Cq, a) — (4, 0.06) 
and conclusion that for sufficiently high cutoff parameter, 
amplitudes are independent of artificial viscosity. 



3.2 Sensitivity to numerical mesh 

We found that amplitudes are almost insensitive to numer- 
ical mesh in case of evolved models, and only slightly sensi- 
tive for models in close vicinity of the blue edge of the insta- 
bility strip. This is illustrated in Tabled First two columns 
give mesh parameters: number of zones, N, and tempera- 
ture of the innermost zone, Tin- In the next columns there 




-5.0 

0.0 0.2 0.4 0.6 0.8 1.0 

phase 



Figure 3. Effects of different artificial viscosity parameters on 
bolometric light curve. Calculations performed for test model C 
(Fig.HJ. 



are bolometric light amplitudes and percentage amplitude 
difference in comparison to model with standard parameters 
(N=120, Tin = 3 x 10 7 K). Stronger sensitivity of amplitudes 
to numerical mesh in direct vicinity of the blue edge is not 
surprising, since near the edge balance between damping and 
driving is a delicate function of our relatively coarse mesh. 
Test models with different location and value of anchor have 
amplitudes varying on similar level as presented in Table |2] 
The chosen depth of the envelope (Tin = 3 x 10 r K), seems 
sufficient to obtain reliable amplitudes. We also note that 
linear properties agree well with result s obtained with mor e 
accurate linear code of Dziembowski l|Dziembowskil IT977f ) . 
We conclude that for most of the main and post-main se- 
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Figure 2. Non-linear limit cycle work per zone for different artificial viscosity cutoff parameters: a = 0.02 (left panel) and a = 0.04 
(right panel). Calculations performed for test model C (Fig.flJ. Surface at right. 



Table 2. Bolomctric pcak-to-peak amplitude, A bo i, sensitivity to 
mesh parameters. N stands for number of zones and T; n corre- 
spond to inner temperature of the model envelope. 



ABC 
N T ln A hol A hol A hol 

120 3 x 10 7 0.50 - 0.95 - 1.05 

160 3 x 10 7 0.56 +12% 0.98 +3.2% 1.06 +1.0% 

120 4 x 10 7 0.42 -16% 0.93 -2.1% 1.03 -1.9% 



always negative, thus only the fundamental and first over- 
tone modes are considered in the following. Assuming that 
model sequences are non-resonant (as will be discussed in 
Section 5, parametric resonance, 2oji ~ wo + W2, play a role 
here, but only in a very narrow range of logT^), and may 
be described by cubic amplitude equations, fo ur scenarios 
are possible (e.g. iDziembowski fc Kovacslll98H ). If 70,1 < 
and 71,0 > then the first overtone (IO) is the only attrac- 
tor of the system. If reverse is true (70,1 > and 71,0 < 0) 
fundamental mode (F) is the only attractor. If both limit 
cycles are stable, 70,1 < and 71,0 < 0, either fundamen- 
tal or first overtone pulsation is possible (so called either-or, 
E/O, situation). Finally if both limit cycles are unstable, 
70,1 > and 71,0 > 0, star must undergo double-mode (DM) 
pulsations. It is clearly visible in Figure [4] that fundamen- 
tal mode pulsation is dominant. First overtone pulsation is 
restricted to the vicinity of the blue edge. Between the fun- 
damental and first overtone pulsation domains, E/O or DM 
(10 Mq) domains emerge. In case of Z = 0.015 only funda- 
mental mode pulsation is possible. In Figure [S] we display 
the F-band peak-to-peak amplitudes as a function of model 
location on the HR diagram. For E/O models and one DM 
model, only the fundamental limit cycle amplitude is plot- 
ted. In Figure[6]amplitudes are plotted vs. the model period. 
Only main sequence models (and one model located directly 
after the TAMS) are included in Figure H3 For E/O models 
amplitudes of both limit cycles are plotted. Open circles 
in Figure [S] correspond to observed monoperiodic variables. 
Their periods and peak-to-peak F-light amplitudes are col- 
lected in Table [1] Several conclusions may be drawn from 
Figures [S] and [5] 

(i) Amplitudes of higher metallicity models (Z — 0.02) 
are larger than those of lower metallicity models (Z — 
0.015). It is a simple consequence of k- mechanism acting 
in the metal opacity bump. 

(ii) The highest amplitudes are obtained for intermediate 
mass models (10-11 Mq), while saturation amplitudes are 
smaller for less massive models (8-9 Mq) and there is a con- 
tinuous fall of amplitudes for higher masses (M > 11 Mq), 
supported by additional calculated sequence of 20 Mq (Z — 
0.02). Near the TAMS and in the post-main sequence phase, 
saturation amplitudes are smaller than in the middle of main 



quence models, limit cycle amplitudes are reliable to within 
3 per cent. 



4 RESULTS 

In this Section we present results of main survey (Sec- 
tion 4.1) and their interpretation (Section 4.2). In Sec- 
tion 4.3 we discuss the effects of different opacities and in 
Section 4.4 justify the neglection of convection. 

4.1 Main survey (OPAL opacities) 

Location of the calculated non-linear models on the theoreti- 
cal HR diagram is shown in Figure|3]for Z — 0.02 (left panel) 
and Z = 0.015 (right panel). Since evolution becomes very 
fast shortly after the TAMS, the TAMS is considered as an 
effective red edge of the pulsation domain. Thus, only a few 
models are calculated beyond the TAMS. Figure [4] presents 
modal selection information, obtained through linear anal- 
ysis and Floquet analysis of the non-linear limit cycles. If 
only one mode is linearly unstable, then this mode is the 
only attractor of the system. If both modes are linearly un- 
stable, selection scenario may be drawn from the Floquet 
analysis of the fundamental and first overtone limit cycles, 
which provides switching rates (non-linear growth rates, Flo- 
quet stability coefficients) toward the first overtone, 71,0, 
and toward the fundamental mode, 70,1, respectively. Gen- 
erally, 7ij are a measure of the linear stability of the limit 
cycle of mode j to perturbation with respect to all the pos- 
sible modes i. In all our calculations 702 1 and 7i>2 are 
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sequence band, but this fall of amplitudes is not very signif- 
icant (Fig.©. 

(iii) The model amplitudes are higher than the observed 
ones (Fig. HJJ . Considering rough averages of all the model 
and observed amplitudes (Table [TJ, the former ones are 
higher approximately by factor ~4 and ~6 in case of Z — 
0.015 and Z — 0.02 model sequences, respectively. Calculat- 
ing average amplitude for non-linear models we used main 
sequence models of masses 9-14 Mq , as these models corre- 
spond to observed variables (cf. Fig.[T]). For E/O models and 
one DM model we chose the fundamental mode amplitudes 
(choosing first overtone amplitudes instead, has no effect on 
presented estimation). 

Calculated single mode saturation amplitudes are signifi- 
cantly higher than amplitudes of monoperiodic f3 Cephei 
variables. As is clearly visible in Figures [5] and [6] theoret- 
ical amplitudes may by lowered be decreasing metal abun- 
dance of the models. After lowering metal abundance from 
Z — 0.02 to Z — 0.015 theoretical amplitudes are still signif- 
icantly higher than the observed ones. At the same time the 
instability strip shrinks and leaves a lot of stars beyond the 
blue edge, as is well visible in Figure Q] Lowering Z further 
will increase this discrepancy. Thus, by lowering the metal 
abundance Z, we are not able to match simultaneously the 
observed amplitudes and the instability domain. 

The natural explanation for described amplitude dis- 
crepancy is, that additional, undetected modes must be ex- 
cited in apparently monoperiodic variables in order to ac- 
count for the low amplitude of the only observed mode. 
Since these additional modes are not visible, either their 
amplitudes are very low or they are of high-degree I, and 
are not visible due to averaging. In the following we show 
that the collective instability saturation by acoustic modes 
is sufficient to explain the amplitudes observed in (3 Cephei 
variables. 
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Figure 6. V-band amplitudes of the non-linear models plotted 
against their period. In case of E/O models both the fundamental 
and first mode amplitudes are plotted. First overtone models (10— 
14 Mq) concentrate near periods of 0.14 days. Circles correspond 
to the observed, apparently monoperiodic /3 Cephei stars. Upper 
panel for Z = 0.02 and lower panel for Z = 0.015. 



4.2 Collective saturation of the pulsation 
instability 

Collective saturation of the pulsation instability may be de- 
scribed by amplitude equations. For each i-th mode of the 
n non-resona ntly interacting modes amplit ude changes are 
given by (e.g. iDziembowski fc Kovacslll984r ): 

^■=ji(l-a a Al + ...-a in Ai)Ai (1) 

where onj are saturation coefficients. Pulsation instability 
is saturated when all dAi/dt — 0. Thus, it is clear, that 
collective saturation of the pulsation instability by several 
acoustic modes will lead to amplitudes lower than in single 
mode case. Quantitative estimates may be drawn assuming 
that all acoustic modes are similar, that is have the same 
saturation coefficients, = a. Then, a simple consequence 
of equations (Jl| is that pulsation instability may be satu- 
rated either by a single mode with amplitude A = \j\fa 
or by n simultaneously excited modes, each of amplitude 
A n — 1/y/na (thus, lower by a factor %/n than in single 
mode case) . In the following, we consider average model and 
average monoperiodic /3 Cephei variable. As was already dis- 
cussed in the previous section, for Z = 0.02 model ampli- 
tudes are ~6 times higher than observed amplitudes (~4 



times for Z = 0.015). Thus with ~36 modes of the same 
amplitude as is observed in average monoperiodic variable, 
one is able to saturate the pulsation instability (with ~16 
modes in case of Z = 0.015 models). Of course these are 
lower limits for number of modes, since in this hypothetic 
multiperiodic variable all the modes should be visible. We 
stress that these are rough estimates, since saturation coeffi- 
cients of different modes do differ. However, it is clear that if 
collective saturation is responsible for amplitude limitation 
a dozen or so modes must be available (linearly unstable). 
In order to match the amplitudes, collective saturation also 
requires higher number of excited modes for more metal 
rich varia bles. Fortunately this is the case as linear the- 
ory shows l|Dziembowski fc Pamvatnvkhll993h . We also note 
that for observed multiperiodic (3 Cephei variables, higher 
numb er of detected modes is correlated with high er metal- 
licity ijDaszvnska-Daszkiewicz fc Niemczurall2005h . 

To check whether a dozen or so non-radial modes 
are available (linea rly unstable) we used LNA code of 
iDziembowskil l| 19771 ) . In this code, full evolutionary stellar 
model is subject to linear pulsation analysis. We determined 
the number of linearly unstable acoustic modes for models 
of different masses (Z = 0.02), located in the centre of the 
main sequence band. We considered modes of degree I < 17 
and dimensionless frequency, a (a = lu/ ^/AirGp) in a range 
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1 < (J < 5. This choice is arbitrary, but well separates groups 
of acoustic, gravity and mixed modes (separation is clear for 
lower masses, and worsens for higher) . Results are presented 
in Table|3j where number of unstable modes, taking into ac- 
count rotational splitting, is given. One should look at orders 
of presented values and trends among them, rather than at 
specific values, since single additional, high degree mode, 
may affect the results presented in the Table. As we have 



checked, the number of unstable modes doesn't vary much 
along evolutionary track. However, the trend of increasing 
number of modes with increasing mass is clearly visible. It 
is clear that considering only acoustic modes, more than 
enough modes are available for collective instability satu- 
ration. We note that not all linearly unstable modes must 
take pa rt in the saturation , as was shown in case of S Scuti 
models |Nowakowskill2005h . The number of linearly unstable 
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M 


n 


8M 


56 


9M Q 


159 


10 M 


278 


11 Mg 


385 


12 M 


418 


14 M 


475 


20 M Q 


665 



Table 3. Results of detailed linear analysis of selected models, n 
stands for the number of unstable modes considering rotational 
splitting. 



modes is also much higher, than observed in multiperiodic 
(3 Cephei stars. Thus, assuming that each number given in 
Table [3] is representative for entire evolutionary track of a 
given mass and that only one third of these modes take part 
in the saturation process, we rescaled theoretical single mode 
amplitudes (Figure [6}, under assumption of collective satu- 
ration. Results are presented in Figure [7] Using only part of 
the linearly unstable acoustic modes, we lowered predicted 
amplitudes to the observed level. Thus, we argue that collec- 
tive instability saturation is sufficient to explain amplitude 
limitation in the f3 Cephei pulsators. 

A closer look at Figure [7J shows that the strongest re- 
duction of amplitudes occurs for 20 Mq track. However, pre- 
dicted amplitudes are still above current detection limit. 
Thus, under our simplifying approximations, non-linear ef- 
fects cannot explain lack of massive and luminous /3 Cephei 
pulsators. Fortunately, three candi dates for such f3 Cephei 
pulsators were found very recently (|Pigulskill2007l) . We note 
in passing that evolved massive pulsators of a Cygni type 
are k nown to occur in th e upper part of the (3 Cephei do- 
main (jLefever et alj|2007f ) . They are located near the TAMS 
and beyond the TA MS. Their variability is interpreted by 
iLefever etT al. (2007) in terms of g-mode pulsation. As such, 
these varaibles are not genuine /3 Cephei pulsators, but 
rather high-mass analogues of the SPB stars. Our calcula- 
tions, which are directly applicable only to p-mode pulsators 
cannot be used to interpret amplitudes of these stars. 

The lack of pulsators at the low luminosity end of the 
theoretical instability strip, can not be explained by non- 
linear effects, either. Even under assumption of collective 
saturation, predicted amplitudes are well above the detec- 
tion limit. However, we stress that these conclusions rely 
on a simplifying assumption about similar properties of all 
acoustic modes. This simplification is useful, since it allows 
for some quantitative conclusions, but in real star satura- 
tion coefficients of modes do differ. One should remember 
that transformation from observational to theoretical HR di- 
agram is uncertain. Possible systematic error in temperature 
transformations may explain lack of less massive pulsators 
near the TAMS, and presence of a few pulsators of higher 
mass beyond the blue edge (see F igure Q. Also using a new 
solar mixture of lAsplund et all (2005) partly solves these 
problems as it leads to theoretical in stability strip more con- 
sisten t with the observational one (jPamvathykh fc Ziomekl 
l2007f) . 

As we have shown, collective instability saturation 
is quite sufficient to explain the observed amplitudes of 
(3 Cephei variables. A possible drawback of this mechanism 
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Figure 7. Theoretical amplitudes, recalculated under assumption 
of collective saturation of the instability mechanism by 1/3 of 
linearly unstable acoustic modes. In addition to monopcriodic 
j3 Cephei variables (open circles) also multiperiodic are plotted 
(filled circles). 



is that excitation of many high-Z modes results in broadening 
of spectral lines, which might be too high when compared 
to observations. Thi s problem has already been discussed by 
iDziembowskl (|l980h in context of 8 Scuti variables. Here we 
follow his reasoning. 

Excitation of a pulsation mode induces time-dependent 
broadening and distortion of spectral lines. The amplitude 
of this effect is given, roughly, by amplitude of the photo- 
spheric velocity variations. In case of collective instability 
saturation, time-dependent contributions from hundreds of 
modes with different periods add up to form a line profile. 
This would mimic macroturbulence and would contribute 
to the measured width of the spectral line. Exact modeling 
of this "collective" line-broadening is a very laborious task 
and is beyond the scope of this paper. Here, we estimate 
this broadening by estimating the root mean square macro- 
turbulence velocity induced by all pulsation modes involved. 
Again, we assume that all modes have the same saturation 
coefficients and that they can be treated as statistically in- 
dependent. Single radial mode creates photospheric velocity 
field vo{t) with amplitude V , which is given by our nonlin- 
ear calculations. V of order of 100 km s _1 is typical for our 
models. When pulsation instability is saturated by n acous- 
tic modes, amplitudes of individual modes are a factor of 
\fn lower, namely v„(t) ~ V / ^/ncosui n t (see previous dis- 
cussion). The root mean square macroturbulence velocity 
can be estimated by simple averaging: 



<5>(*)> 



> < COS^ LJ„t > 



V_ 

71 



(2) 



Thus, on the basis of our nonlinear calculations, we ex- 
pect the line broadening in collective saturation scenario to 
be of the order of 100 km s _1 . This number is, generally, 
comparable to v sin i value s measured in (3 Cephei variables 
l|Stankov fc Handler! [20051 . fig. 2), however, for many stars 
the observed line broadening is smaller. This poses a poten- 
tial problem to the proposed scenario. A possible way out 
of this difficulty is to include the g-modes (omitted in our 
analysis) into the saturation process. /3 Cephei model cal- 
culations kindly provided by Dziembowski (private comm.) 
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show, that photospheric velocity variations of the g-modes 
should be smaller than those of the p-modes. This is true, 
assuming the same flux perturbation amplitudes in the driv- 
ing zone, which is the natural normalisation when instability 
saturation is discussed. 



4.3 Effects of different opacities 

Linear properties of the (3 Cephei models are sensitive to 
opacities being used. Particularly, the OP o pacities lead 
to wid er instability strip, as demonstrated by IPamvathvkhl 
(1999) (during preparation of this paper, we have noticed 
that also different methods of interpolation of the OPAL 
opacities, lead to slightly different results, but since this ef- 
fect is small, and has no effect on presented conclusions, 
we won't discuss it). We calculated one sequence of models 
(M = 11 M Q , Z = 0.02) using the OP opacities. Models 
have the same log L and log T e g parameters as models us- 
ing the OPAL opacities, since we are interested in opacity 
effect on pulsation properties. Models with modal selection 
information are presented in Figure [5] To allow comparison, 
evolutionary track with the OP results is artificially shifted. 
Vertical ticks mark the edges of the fundamental mode in- 
stability strip. One may observe completely different modal 
selection scenario, with wide region of DM models for the 
OP sequence. Since robust DM behaviour, in purely radia- 
tive models, was not en countered up to date (cf. however 
iKovacs fc Buchler|[l9 93) , we will discuss it in Section 5, here 
concentrating on amplitudes only. Figure [9] presents peak- 
to-peak F-band amplitudes of the fundamental limit cycle 
as a function of mo del's effective temp erature. Size of the 
circles correspond to IStellingwerj (|l97&t ) linear growth rate, 
rj. It characterises the envelope driving: with rj = 1 the en- 
tire envelope drives, while with r\ — — 1 the entire envelope is 
dissipative. Amplitudes obtained with the OP opacities, are 
generally higher, which is connected with higher rj values. 
However for temperatures higher than log Toff « 4.37 visible 
difference between the OP and OPAL results is connected 
with the location of the blue edge, significantly shifted to 
the blue in case of the OP opacities. For lower temperatures 
amplitudes obtained with the OP opacities are compara- 
ble to amplitudes obtained using the OPAL opacities. Thus, 
presented conclusions remain unchanged. 



4.4 Effects of convection 

Contrary to the classical pulsators, convection is not ex- 
pected to play significant role in hot f3 Cephei stars. As 
we checked for 11 Mq, evolutionary track calculated with 
mixing-length convection (a = 1) overlap the track calcu- 
lated without convection (see also IPamvathvkhl Il999). Im- 
pact of convection on pulsation models was checked us- 
ing direct time integration convective code. We used La- 
grangea n code with time-depend ent convection model of 
iKuhfud ( 1986T I in the version of IWuchterl fc Feuchtingerl 
(199J). All other hydrodynamics and procedures are the 



same as were used in radiative hydrocodes. Since growth 
rates of j3 Cephei stars are small, convergence to limit cycle 
may require calculation of more than 10 5 pulsation cycles. 
For this reason only four models and two sets of convective 
parameters were checked (a G {1.0, 0.5}, a m = 0.4; a s , cd, 
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Figure 8. 11 Mq model sequence calculated using both the 
OPAL and OP opacities in the pulsation codes. Since evolution- 
ary track was calculated using the OPAL opacities, track with 
the OP results is artificially shifted to allow comparison. Vertical 
ticks mark the linear edges of the fundamental mode instability 
strip. Arrows indicate the position of the 2u>i ~ u>o+u!2 resonance 
centre. 
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tudes of 11 Mq models calculated with both the OPAL (solid line) 
and OP (dashed line) opacities. Size of the circles correspond to 
model's Stellingwerf 's rj. 



-y r - standard values as given by IWuchterl fc Feuchtingerl 
1998). Turbulent pressure and overshooting were neglected 
(a p = at = 0). Tested models are evolved and massive 
models pulsating in the fundamental mode. These models 
are characterised by relatively high growth rates (limiting 
amplitude is achieved faster) and have stable first overtone 
(transient phase is short). Models are marked by crosses in 
Figure [1] As initial model for integration we used fully con- 
sistent static convective envelope, constructed in the same 
way (number of zones, anchor, depth of the envelope) as for 
radiative models. We find that during the contraction phase 
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M 


^bol 










a = 1 


a = 0.5 


14 M 


0.85 


0.83 


0.85 


12 M Q 


0.98 


0.95 


0.98 


11 Mq 


1.05 


0.99 


1.05 


10 M Q 


1.16 


0.98 


1.14 



Table 4. Comparison of bolometric light peak-to-peak ampli- 
tudes of selected radiative j) and convective (A£ j) mod- 
els computed with mixing-length equal 1 (third column) and 0.5 
(fourth column). 

up to 10 per cent of the flux may be carried by metal-bump 
induced convective zone. Limiting pulsation amplitudes are 
lower than in case of radiative models, but the difference 
is generally small, which is shown in Table [4] Lower ampli- 
tudes for convective models are expected, since part of the 
flux is carried by convection and thus the K-mechanism is 
less effective. We note up to 18 per cent difference between 
radiative and a = 1 convective amplitudes for less massive 
models. With a = 0.5 convective and radiative amplitudes 
are almost the same. We argue that a = 0.5 is a better value, 
since the mixing-length should be smaller than the extent of 
the convective zone. This condition is satisfied better with 
smaller a values. However, we stress that even with a = 0.5, 
at some pulsation phases the extent of the convective zone 
is comparable with the mixing-length. Thus, this model of 
convection, working well in case of classical pulsators, is at 
the limits of its applicability here and is used only for estima- 
tion of the impact of neglecting convective energy transfer 
on calculated amplitudes. We also expect that convective 
effects get weaker as one moves from the TAMS (where our 
test models are located) toward the ZAMS. Construction of 
the static convective envelopes for 10 Mq evolutionary track 
clearly shows this. For this sequence maximal static convec- 
tive flux for model presented in Tablets 2.7 per cent of the 
total flux, while for the model closest to the ZAMS it is less 
than 0.4 per cent. Of course dynamical effects lead to higher 
values of convective flux in non-linear models, but we expect 
the same trend. We conclude, that although amplitudes of 
convective models are systematically lower, this difference is 
not significant for our discussion. 



5 DOUBLE-MODE MODELS 

Although DM RR Lyrae and Cepheids are quite com- 
mon, for many years purely radiative codes failed to re- 
produce the DM b ehavi our. In case of RR Lyrae models 
iKovacs fc BuchieU (|l993l ) showed, that for particular range 
of artificial viscosity parameters, the DM state is possible, 
but still their models were strongly dependent on numerical 
details, and were not in satisfactory agreement with obser- 
vations. It was the inclusion of convecti on, that finally led 
to robust and sati sfactory DM models (|Feuchtingerl 1 19981 ; 
iKollath et al.ll 199*81 ). To our surprise we encountered one DM 
model in our radiative models built with the OPAL opaci- 
ties (Fig. [4| , and a wide region of DM models in radiative 
models that use the OP opacities (Fig. |8}. This was un- 
expected, since among observed multimode /3 Cephei stars, 
pulsat ors with two radial mod es only have not been detected 
so far (|Aerts fc De C at 2003). Linear analysis showed, that 



despite obvious role of different opacities, also a paramet- 
ric resonance, 2lji ~ luq + u>2 , may be connected with DM 
phenomena. The resonance centre is indicated by the long- 
dashed line in Figure [4] and by arrows in Figure [8] As will 
be shown in the following, different mechanisms are respon- 
sible for DM-behaviour in case of the OPAL and the OP 
opacities and thus, models computed using both opacities 
are discussed in separate subsections. We also stress that 
our modal selection analysis relies only on Floquet stabil- 
ity coefficients (switching rates). If stability coefficients of 
both the fundamental and the first overtone limit cycles are 
positive, neither limit cycle is stable, and DM behaviour is 
unavoidable. However, only direct time integration provides 
information about amplitude stability of such solution. Due 
to very small growth rates, integration covering several hun- 
dred thousand cycles would be required. Such a calculation 
is prohibitively time-consuming and was not performed here. 

5.1 The case of the OPAL opacities 

It is clearly visible in Figure [4] that the only DM model 
found in the main survey (10 Mq), is quite far from the res- 
onance centre and thus, is not connected with it. It also 
seems, that the resonance does not lead to any change in 
the modal selection. However, this is not the case. We have 
calculated additional models of 10 Mq in direct proximity 
of the resonance centre, and close to the only DM model 
visible in Figure [4] Stability results are presented in Fig- 
ure [TO] Growth rates (both linear and non-linear) are plot- 
ted against parameter A, defined as A = 2wi/(wo + ljz). A 
is calculated using linear frequencies. The shaded bar indi- 
cates modal selection. Two DM domains emerge, with two 
different mechanisms being responsible: 

(i) Resonant DM. The most striking feature visible in 
Figure [TU] is the prominent peak of switching rate toward 
the fundamental mode, 70,1, with no doubt connected with 
the resonance. Otherwise negative, 70,1 becomes positive at 
the resonance centre. Thus, the discussed resonance leads to 
the destabilization of the overtone limit cycle. As a result, 
the DM island appears in the middle of the first overtone 
pulsation domain. The domain of the resonant DM models 
is very narrow (AT c fr ~ 200K), and thus, no such models 
are visible in Figure [4] since they fall in between the consec- 
utive first overtone models. We also checked that for 11 Mq 
evolutionary track there are DM models in direct proximity 
of the resonance, very close to the linear fundamental mode 
blue edge. For the 12 Mq evolutionary track, resonance cen- 
tre falls in between the linear edges of the fundamental and 
first overtone modes and the overtone limit cycle is stable 
(70,1 < 0). Thus, the resonant DM models are restricted to 
intermediate masses, since for masses slightly above 12 Mq 
resonance centre falls beyond the instability strip, while for 
smaller masses, below ~9Mq, 71,0 < 0, and thus the funda- 
mental mode is the only stable attractor (Fig. [4]). 

(ii) Non-resonant DM. The second DM domain, at 
A « 1.01 is not connected with any resonance. DM domain 
separates the first overtone and fundamental mode pulsation 
domains. Such behaviour is expecte d in case of non-resonant 
coup ling of pulsation modes (e.g. Dzicmbowsk i fc Kovacsl 
11984) . The non-resonant DM behaviour is also restricted 
to intermediate masses (~10Mq), since it requires that the 
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switching rates do intersect, and are positive at the inter- 
section. For less massive models the switching rates do not 
intersect, while, for massive models, intersection occurs at 
the negative value of switching rates, which leads to E/O 
behaviour. 



5.2 The case of the OP opacities 

In case of the OP opacities, the domain of the DM-behaviour 
is much wider than in case of the OPAL opacities, as is well 
visible in Figure[8] It is also visible that the DM models seem 
to be concentrated near the resonance centre. However we 
will show that the resonance does not play any role here. To 
allow direct comparison with the OPAL results (Figure [10)) 
we calculated several models of 10 Mq in the proximity of 
the resonance. Switching rates are presented in Figure [11] 
(qualitatively the same picture is true for 11 Mq evolution- 
ary track for which modal selection results are presented in 
Figure[8]). Both linear growth rates are positive (the red edge 
of the first overtone pulsation falls at A w 1.018), and for 
clarity are not shown in the Figure. The prominent peak of 
70,1 is well visible, but the resonance producing this peak is 
not responsible for the destabilization of the overtone limit 
cycle since, opposite to the OPAL case, 70,1 > almost ev- 
erywhere. Thus, it is the switching rate toward the first over- 
tone, 71,0 (i.e. stability of the fundamental limit cycle) that 
determines the edges of the DM domain. 71,0 seems to peak 
near the resonance centre and there is a secondary maximum 
at around A « 1.008. This second DM domain is probably 
the numerical artefact (probably connected with the strong 
shock propagation). We argue that the first DM domain falls 
at the A ~ 1 by accident, and thus 2lui ~ ujo +LU2 resonance 
is not responsible for the destabilization of the fundamental 
limit cycle. Our arguments are as follows: 



(i) iKovacs fc Buchlerl (|l993) analysed the 2o»i ~ u>o + ^2 
resonance, using lowest non-trivial order of the amplitude 
equations. They showed that the resonance affects only the 
stability of the overtone limit cycle, 70,1, while it doesn't 
affects the stability of the fundamental limit cycle, 71,0- 

(ii) The same may be shown by numerical experiment. 
To check the effect of the resonance on stability coefficients, 
we calculated additional sequence of 10 Mq , with shallower 
envelope, inner temperature being, Tin = 2.4 x 10 7 K. Re- 
sulting mode frequencies are slightly different in comparison 
to models with standard envelope and thus, the resonance 
centre is shifted to higher effective temperatures. Stability 
results are presented in Figure 1121 The prominent peak of 
70,1 is still well visible at the resonance centre. However the 
maximum of 71,0 has nothing to do with the resonance. 

We note that even if one treats the second DM domain 
visible in Figure [TT] as a numerical artefact, the remain- 
ing picture is far from the one expected in case of non- 
resonant mode coupling, where DM domain separates the 
first overtone and fundament al mode pulsation domains 
ijDziembowski fc Kov acs 1984). Here, the double- mode do- 
main appears as an island in between the fundamental mode 
pulsation domains. We are not able to give a specific reason 
for the destabilization of the fundamental mode visible in 
Figure [TT1 
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Figure 10. Linear growth rates (thin lines) and non-linear 
growth rates (switching rates, thick lines) plotted against A, 
A = 2a>i/(wo +^2)1 for 10 Mq evolutionary track computed us- 
ing the OPAL opacities in pulsation calculations. The resonance 
centre (A = 1) is located at logT off ss 4.378. 
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Figure 11. Non-linear growth rates (switching rates, thick lines) 
plotted against A, A = 2uii/(ioo +^2), for 10 Mq evolutionary 
track computed using the OP opacities in pulsation calculations. 
Arrow marks the location of the ZAMS. Both linear growth rates 
are positive almost in the entire range of A (both growth rates 
are positive at the ZAMS and the red edge of the first overtone 
pulsation falls at A as 1.018) and therefore are not plotted. The 
resonance centre (A = 1) is located at logT e ff fa 4.379. 
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Figure 12. Non-linear growth rates (switching rates) plotted 
against A, A = 2wi/(wo +1^2), for 10 Mq evolutionary track 
computed using the OP opacities in pulsation calculations and 
shallow envelope model. Arrow marks the location of the ZAMS. 
Resonance centre is located at logT c fj = 4.392. 

6 CONCLUSIONS 

In the present analysis we have focused on amplitude sat- 
uration in P Cephei models. We calculated several evolu- 
tionary tracks of two different metallicities, Z — 0.02 and 
Z — 0.015, and masses covering the observational insta- 
bility domain. More than hundred models were subject to 
non-linear pulsation analysis. Radial modes, treated as rep- 
resentative for al l acoustic os c illatio ns, were converged to 
limit cycle using IStellingwerJ l|l974h relaxation technique. 
We have shown that our models are robust, independent 
of numerical treatment. Particularly, contrary to the clas- 
sical pulsators models, for physically correct choice of ar- 
tificial viscosity parameters the saturation amplitudes and 
light-curve shapes of our models are independent of artificial 
viscosity. 

Examination of saturation amplitudes along the model 
sequences, as well as comparison of these amplitudes with 
amplitudes observed in /3 Cephei stars leads to several 
conclusions. Most important concern amplitude limitation 
mechanism in f3 Cephei stars. Single mode saturation am- 
plitudes are significantly higher than amplitudes observed 
for monoperiodic /3 Cephei variables. Assuming similarity 
of all the acoustic modes, we showed that collective insta- 
bility saturation easily lowers predicted amplitudes to the 
observed level. The straightforward conclusion is, that in 
fact all /3 Cephei pulsators should be multiperiodic. On- 
going decrease of detection limit should reveal additional 
modes in monoperiodic /3 Cephei variables. This is fully sup- 
ported by recent progress in intensive multisite campaigns 
leading to the detect i on of more a nd more pulsation modes 
l|Handler et al l |2004 120051 . l2006h . A possible weakness of 
the collective saturation scenario is that the estimated line- 
broadening induced by excitation of hundreds of pulsation 



modes might be higher than observed in some stars. We ar- 
gue that this difficulty can be solved by allowing g-modes to 
take part in the saturation. We also addressed the problem 
of non-uniform filling of the theoretical instability strip, spe- 
cially the lack of massive and luminous p-mode pulsators. 
Although our calculations show decrease of saturation am- 
plitudes while moving to higher masses, we found this effect 
too weak to explain scarcity of such pulsators. 

We have found that our results depend somewhat on 
the opacities being used (OPAL or OP). Differences are sig- 
nificant but do not change our conclusions. A few non-linear 
models were c alculated with turbulent convection based on 
iKuhfufj (|l986h theory. In agreement with expectations, we 
found that convection may be neglected in /3 Cephei models. 

In several of our radiative models we have found a 
double-mode behaviour, with radial fundamental and first 
overtone modes simultaneously excited. This form of pul- 
sation is encountered only in intermediate mass models 
(10-11 Mq). It is numerically very robust and occurs both 
with OPAL and with OP opacities. Depending on the 
specific model, the origin of double mode pulsation can 
be traced to one of two different mechanisms: either to 
the nonresonant coupli n g of the two excited modes (e.g. 
iDziembowski fc Kovacsl (jl98J)), or to the 2lj± ~ uq + LU2 
parametric resonance. 

Admitedly, our radial double-mode models have no di- 
rect application to real stars, because no purely radial multi- 
mode ft Cephei pulsators are known. Nevertheless, our mod- 
els demonstrate how easily low-order parametric resonances 
can destabilize a single-mode radial pulsation. When all lin- 
early unstable (or marginally stable) non-radial modes are 
taken into account, many more resonances of this type can 
be found. This might explain why high ampliude monoperi- 
odic radial pulsators (similar to BW Vul) are so rare among 
/3 Cephei stars. 
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